Properties of high-Tc copper oxides from the nearly-free electron model. 
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The generic band structure of high-Tc copper oxides is simulated by the nearly free-electron 
model (NFE) in two dimensions (2-D) with parameters from band calculations. Interaction between 
phonons and spin waves will cause potential modulations and pseudogaps, and the strength of the 
modulations, the wave lengths and the doping, are all related. A Fermi-surface "arc" is found for 
dynamic spin/phonon waves. The confinement of superconductivity between two limiting dopings 
can be a result of competition with the pseudogap at low doping and weak coupling at high doping. 
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An important structural feature of high-Tc materials 
is the stacking of almost 2-dimensional (2-D) CuO planes. 
All high-Tc cuprates have at least one of these layers in 
the unit cell, and they make the band structure fairly 
simple with one or more M centered Fermi surface (FS) 
cylinders. Band calculations and photoemission agree es- 
sentially on this fact, except for details and for undoped 
materials, which often are antiferromagnetic (AFM) in- 
sulators [lj]. The density of states (DOS) at the Fermi 
energy (Ep) is mainly of Cu-d character with some O-p 
admixture. The views on the mechanism for supercon- 
ductivity and normal state properties diverge. Extensive 
experimental works show that not only the high Tc, but 
also many normal state properties, are unusual. These 
include pseudogaps and magnetic fluctuations, both with 
rich dependencies as function of doping, x, and tempera- 
ture, T. The importance of phonons is evidenced by iso- 
tope effects on Tc and pseudogaps @, Q , and magnetic 
fluctuations are detected by neutron scattering 0]. 

Here is presented a NFE model for high-Tc proper- 
ties with use of parameters coming from ab-initio band 
calculations for pure and hole doped high-Tc materials 
0, Q . These band calculations show important coupling 
between phonon and spin waves. The difficulty with 
ab-initio calculations is that a unit cell with interesting 
phonon displacements and/or spin waves needs to be very 
large. In addition, the bands are not obtained within the 
normal Brillouin Zone (BZ), but within the down-folded 
zone. Such calculations are so far limited to spin waves 
and phonons along one direction only, while in reality 
one can expect modulations along x and y (checker board 
modulations rather than stripes). The qualitative results 
from the one-dimensional (1-D) NFE model and ab-initio 
calculations are the same, and it is worthwhile to extend 
the NFE-model to 2-D. 

A periodic potential perturbation, V(x) = 
Vc.exp(—iQ ■ x), will open a gap of size 2Vq at the 
zone boundary, k = Q/2, in the 1-D NFE bands @, 0. 
This well-known NFE result can be applied to the AFM 
spin arrangement on neighboring Cu along [1,0,0] (with 
wave vector Q) in undoped insulating cuprates. An 
additional modulation, with wave vector q < Q, modifies 
the potential, V(x) — Vc.exp(—iQ ■ x)exp(iq • x). The 
gap moves away from Q/2 to (Q — q)/2, as for ab-initio 
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FIG. 1: Calculated 2-D NFE Fermi surface for a strong dy- 
namic spin-phonon fluctuation with q— (0.85,0.90). 



bands in supercells of different lengths. The combined 
potential modulation leads to stripe like phonon or 
spin-wave patterns, with short wave lengths (periodicity, 
with wave vector q) at high doping and long ones when 
x — > 0, as shown in previous band calculations Q. 

Simultaneous modulations along x and y are yet too 
difficult for an ab-initio band approach, at least for real- 
istic wave lengths. An extension of the NFE model for 
potential modulations along x and y leads to a 3x3 eigen- 
value problem of the form E — k 2 — k y , E—(k x — Q x ) 2 ~k y , 
and E — k 2 . — (k y — Q y ) 2 in the diagonal, and Vq as 
non-diagonal terms. The G-vectors 0, Q x = Q — q x and 
Qy = Q — Qy ar e considered in the basis with q x and q y 
along x and y. The bands are represented in the reduced 
zone as for the normal unit cell. The band width from 
the T-point up to kp is about 0.125 Ry when VQ,q x ,q y 
and x are all zero and the effective mass is 2.5. The po- 
sition of the van Hove singularity (when the bands touch 
X and Y) corresponds to x « 0.2. Potential modulations 
will open gaps near X and Y, while not much happens 
between T and M Q. 

Superconductivity and pseudogaps appear typically at 
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100K, and we will estimate the VQ-parameters for this 
temperature. The mean value of atomic displacements 
u = yJlZksT '/ K), where the force constant (K) for O 
vibrations is of the order 25 eV/A 2 for thu w 50 meV, is 
close to O.Olao, where ao is the lattice constant. Zero- 
point motion is just a bit larger, but it is not selective to 
one phonon [9(. The mean value of magnetic moment 
fluctuations, mo, is obtained through fcsT = ^Fm 2 , 
where F = d 2 E/dm 2 and E is the total energy. For a 
short spin-phonon wave in HgBa2Cu04 3, Vq is about 
6 mRy for tuq ~ 0.09/is per Cu. 

The origin of Vq is two-fold. Structural distortions of 
phononic origin contribute equally to the two spin com- 
ponents of V(x), while for a spin wave there is an op- 
posite phase of the two spin potentials. A static spin- 
polarized modulation of the potential determines AFM 
order in the undoped case. Phonons, and probably also 
long-wave spin- modulations, are dynamic. From the self- 
consistent band calculations for 'half breathing' phonons 
along [1,0,0] it is estimated that Vq 1 varies from 5 to 11 
mRy when the wave length varies from 4ao to 16ao- This 
is when the atomic displacement is of the order O.Olao, as 
when T » 100K Q. The phonon amplifies a spin wave, 
which has twice a long wave length as the phonon, but 
the two types of waves enforce a common gap at the same 
energy \v$. 

The self-consistent convergence of long- wave spin con- 
figurations is very slow in band calculations. In or- 
der to estimate the spin-polarized part Vq we rely on 
the band results for a short spin wave and do NFE- 
scaling for longer waves. Band calculations on undoped 
Hg2Ba4Cu2C>8 show that an AFM, zero-gap state can 
be stabilized with a spin-splitting of 23 mRy on Cu Q. 
(Larger gaps are possible when using a different density 
functional in band calculations for La2Cu04 [HI). This 
determines the Vq for the basic AFM spin arrangement 
between nearest Cu neighbors. A well-known feature of 



TABLE I: Spin density coefficients (A cos 2 (7r:r/2) and B 
sin 2 (7ra;/2)) for the NFE state below the gap at the zone 
boundary along q x . The exchange splitting for spin fluctu- 
ations, Vq , is estimated from the scaling procedure as de- 
scribed in the text. The last column shows the potential pa- 
rameter for phonons, Vq , which are interpolated from band 
results for Q < 0.9. A saturation is assumed for Q > 0.9. 
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FIG. 2: Calculated DOS for z=0.02 (thin line), x=0.16 (dot- 
ted) and x=0.26 (heavy line). The latter case is for non- 
equivalent modulations along x and y, and a secondary dip 
below Ef is seen. 



the 1-D NFE model is the spacial separation of the states 
below the gap (^ 2 _(x) = sin 2 (Gx/2)) and above the gap 
(*+(£) = cos 2 (Gx/2)) :7ij. Let the state below (above) 
the gap coincide with the Cu with the attractive (repul- 
sive) potential for one spin. The spin density of the first 
(second) spin is given by the sin 2 - (cos 2 )-term, and the 
phase on the nearest Cu neighbor differ by tt. The densi- 
ties feed the exchange splitting of the potentials. If they 
diminish it leads to a smaller exchange, which leads to 
smaller densities and so on, in a self-consistent manner. 
Thus, the spin density is a driving force behind the near 
neighbor AFM configuration, and it depends on the op- 
timal occupation of the two states. 

An additional modulation will reduce the spacial sep- 
aration of the two states at G/2. Table I shows how the 
lower state is mixed as function of Q x . As Vq depends 
on the near-neighbor interaction one can expect a reduc- 
tion of Vq by the coefficients given in Table I. This does 
not include self-consistent feedback, or the closing of the 
pseudogap at T* (see later) , effects which both should de- 
crease V^. Usually T* is well above the 100K at which 
the parameters are estimated. But when T* goes to zero 
at large doping, it will reduce Vq more (the two last 
lines in Table I include rather arbitrary reduction factors 
of 0.9 and 0.8). On the other hand, the coefficients Vq 
should increase because of coupling to phonons for all q > 
0. This effect is estimated to increase the moments by at 
least 30 percent for a short phonon wave with appropriate 
value of u [14j , and the values in Table I include a factor 
of 1.3. It can be recalled that Vq at q x — q y = 0.5 from 
the simple rescaling, 6 mRy, agree with the independent 
estimate from the best converged band calculation for 
that spin wave (4a ) in doped HgBa 2 Cu04 @. 

An example of the FS is shown in figure 1, which is a 
sum of several calculations with Vq ranging from 10 to 30 
mRy in order to simulate dynamical waves. The position 
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FIG. 3: The relation between doping x and q-vectors of the 
modulation. The q-vectors along x and y axe equal for doping 
below ~ 0.18. The modulation along one direction remains 
fixed (8 ao for the spin part) for larger doping, with a weaker 
modulation of longer periodicity in the perpendicular direc- 
tion. 
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FIG. 4: The DOS at E F and T* as function of doping x. The 
DOS is at the minimum of the pseudogap and T* is l/4th of 
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of Ep is chosen at the minimum of the sum of the partial 
DOS functions, and the plot is finally symmetrized along 
the T-M line. The FS remains as a sharp "arc" on the 
diagonal direction, while it is washed out near the limits 
of the zone because of the fluctuations of Vq . The result 
is compatible with the observation of a small section of a 
FS-arc at low T 
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The arc widens for larger T, when 
most of the potential modulations and the pseudogap are 
gone. Static potential modulations will bend the outer 
sections of the FS towards the X, Y — M lines [fj , which 



in a repeated zone looks like a second " ghostband" . 

Calculations of the DOS (fig. 2) and the relation 
between doping and q-vectors (fig. 3) are made using 
Vq = + from Table I as input. The Q-vectors 
are allowed to vary from 1.0 to 0.75 in x and y directions, 
and the doping is optimal when Ep coincides with the 
energy at the DOS minimum, at the pseudogap. Three 
examples of doping are shown in figure 2. At low dop- 
ing (x < 0.18) there is a nearly linear relation between x 
and q (see fig. 3), and the Q-vectors along x and y have 
equal lengths down to the minimum at 0.75 of the zone 
boundary limit. There is no possiblity to obtain a mini- 
mum in the DOS for larger doping with equal | Q x | and 
| Q y |, the values of Vq are too small. However, if one of 
the Q-vectors remains fixed at the value 0.75 (assumed 
to correspond to the shortest possible magnetic modula- 
tion [4j) while the other one increases, it is possible to 
follow the pseudogap further towards large x, with Vq 
coming from Table I. (The larger of the two Q-vectors 
makes a dip in the DOS, but it is weak and below Ep, 
see fig. 2, so its modulation should be harder to detect.) 
The result in Fig. 3 is qualitatively similar to the dop- 
ing dependence observed in La(2- a: )Sr a: Cu04 by Yamada 
et. al. [l3j], although saturation of the periodicity oc- 
curs near x ~ 0.18 compared to about 0.12 as observed. 
Other combinations of Q-vectors and Vq can give larger 
x, but with weaker gaps or unrealistic values of Vq. 

The unperturbed NFE band is perfectly isotropic, but 
the real band structure may have different dispersion in 
different directions. A FS which extends more towards 
the diagonal than towards X and Y, can be modelled by 
an anisotropic effective mass. Here, through multiplica- 
tion of the mass with ((| k x \ + \ k y |)/(| k x + k y |))a 
we test a ~ 10 percent anisotropy. The result is that the 
scale of the doping in figs. 3 and 4 becomes compressed, 
and the breaking point in fig. 3 is at a; =0.13 instead 
of at 0.18, which would fit better to the results by Ya- 
mada et. al. jl3j ]. The values of N(Ep) will go down 
slightly, especially for small x. This shows that details of 
the real band structure can be important for the quanti- 
tative results. Ab-initio bands in undoped HgBa2Cu04 
Q suggest that the FS retracts towards the diagonal. 

The relation between doping and wavelength from the 
(1-D) band calculations for long supercells appears to be 
very exact. For instance, a doubling of a cell will result 
precisely to a factor 2 in the doping. This is because 
practically all Ej~ (near Ep) in the compressed 1-D BZ 
become gapped. Many states in the interior of the zone in 
the 2-D NFE model are not concerned by gaps. All states 
contribute to the DOS and the energies of local gaps do 
not correspond to the pseudogap in the total DOS. One 
could return to the more 'exact' relation if the 2-D sheet 
of the NFE band behaved very rigidly with similar gaps 
everywhere. 

The T-dependent Fermi-Dirac occupation leads to a 
quenching of the spin wave by the feedback of the spin 
density on to the potential. The quenching temperature 
defines T* , which in a 1-D model is about jA/ks, where 
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A w Vq S is the gap at T=0 0. Figure 4 shows T* and 
N(Ep) as function of the doping from the model. 

The change of q as function of coupling strength was 
not forseen in the 1-D band calculations. If this fact is 
true, then one can expect longer wave length for weaker 
coupling, near T* for instance. This implies that the 
gap moves away from Ep , which is less favorable for the 
stability of the wave, and the process of quenching at T* 
may accelerate. The same mechanism predicts isotope 
shifts on the q-vectors of magnetic fluctuations, since a 
heavier mass is expected to decrease Vq through smaller 
u 0. 

The 1-D band calculations show that a phonon and a 
spin wave, which along [1,0,0] differ by a factor of two 
in wave length, tend to open a gap at the same energy. 
However, the spacial shape of the potential perturbation 
caused by a phonon and a spin wave is different, (this is 
less clear from the NFE-model) and regarding supercon- 
ductivity it is not expected that the Vq 's of the two waves 
work together. The spin wave, and equal-spin pairing 
(ESP), is probably most important for superconductiv- 
ity, since Vq can become very large, as when a favorable 
phonon displacement is assisting. These arguments sug- 
gest ESP as a mechanism for superconductivity, but it 
has to compete with the pseudogap which removes states 
and DOS near Ep, more so at low doping, see figure 4. 
This scenario is corroborated by recent femtosecond spec- 
troscopic measurements on cuprates showing competing 
order from something like a pseudogap within the super- 
conducting gap [l5j . The coupling for spin-fluctuations 
decreases towards the over-doped side, as is reflected by 
the decreasing V^-values in Table I. Also the coupling 
to phonons disappears at too short wave lengths, since 



no spin wave shorter than 4ao can co-exist with a "half- 
breathing" phonon. The coupling parameter A ~ N ■ V 2 
and Tq will therefore vanish at the extreme dopings, as 
can be deduced from figure 4 by the low DOS for x — > 
and the small T* (which is proportional to Vq) at large 
doping. The present results are not sufficiently complete 
for an evaluation of Tq through a BCS-likc formula, but 
qualitatively it is expected that the limits of a " X^-dome" 
are shaped by the lines for T* and N(Ep), as in figure 
4. In order to increase Tc on the under-doped side one 
should increase N(Ep). One possiblity is to make the 
2-D sheet of the band less rigid in order to restrict the 
gaps to the exterior of the BZ, perhaps by pressure. 

Many typical high-Tc features, such as FS-arcs, pseu- 
dogaps, T*, and the g(x)-dependence, can be described 
qualitatively by the NFE simulations. The total Vq de- 
termines the size of the pseudogap, but Vq disappears 
above T* and the same gap cannot be supported by Vq 1 

alone. The smaller Vq 1 makes a weak dip in the DOS at 
lower energy, away from Ep and away from optimal dop- 
ing. It is suggested that superconductivity is caused by 
ESP, and Vq, which leads to confinement of Tc between 
two limiting dopings. The parameters are based on pre- 
vious band calculations for phonons and spin waves in 
doped systems, and the values in Table I are the most 
probable Vq to use in the 2D-NFE model. Still, it is 
unavoidable that some of the estimations are very un- 
certain. But finally, the extreme simplicity of 2-D NFE 
model makes it a toy model, where other solutions of the 
parameter space can be tested. 

I am grateful to B. Barbicllini and C. Berthod for var- 
ious discussions. 
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